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ABSTRACT 

We investigate the recently snggested scheme of independent mass matrices for neutrinos 
and antineutrinos. Such a CPT violating scheme is able to account for all neutrino data with 
the three known flavors. For atmospheric neutrinos this means that it is possible to have 
different mass squared differences driving the oscillation for neutrinos and antineutrinos. We 
analyze the atmospheric and K2K data within the simplest scheme of two neutrino oscillation, 
neglecting electron neutrino oscillation. We hnd that the preferred region is close to the CPT 
conserving mass spectra. However the spectra with the antineutrino mass squared difference 
about or larger than 0.1 eV^ and the neutrino mass squared difference about 2 x 10“^ eV^ is 
not signihcantly disfavored. In this parameter region the atmospheric data are independent of 
the antineutrino mass squared difference. Therefore no useful constraint can be put on CPT 
violation effects contributing to different masses for the neutrinos and antineutrinos. 
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1 Introduction 

Many elementary particles, like the electron and the kaons, provide tight bounds on possible 
CPT violating effects contributing to different masses for the particle and its antiparticle. For 
instance for the electron and the positron we have 

|mp+ — Trip- 1 Q ^ P.V 

^^ < 8 X 10"® , CL = 90% . (1) 

^average 

As is well known, CPT conservation implies the equality of the neutrino and antineutrino 
survival probabilities in vacuum [Q, though matter effects can produce fake CPT violating 
effects [^. The atmospheric neutrino data involve both the particle and the antiparticle 
channels and are therefore suitable for a study of possible CPT violation in the neutrino sector. 
The idea to use neutrino oscillation to search for CPT violation was hrst proposed in Ref. [Q. 
Naturally as the atmospheric neutrino experiments are probing mass squared differences and 
not the absolute neutrino mass, these will be the quantities which might be restricted by the 
data. The interest in CPT violation arises due to a recently suggested scheme which is capable 
of solving all neutrino anomalies without the use of a light sterile neutrino |^, 

At present three neutrino anomalies (atmospheric [^, solar p and LSND 0) exist, all 
requiring different Am^’s when interpreted in terms of neutrino oscillation. Therefore a CPT 
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Figure 1: Schematic view of the masses of neutrinos and antineutrinos 


conserving three neutrino framework cannot account for all anomalies. This has also been 
explicitely shown in theoretical hts of the atmospheric data BP- Consequently one has to 
go beyond standard explanations to solve all anomalies. 

A possible solution could be the existence of a light sterile neutrino. Several studies of 
such four neutrino models have been performed and the current situation has been presented 
in Ref. |l^ . The four neutrino models give an acceptable £t when fitting all available data. 


However, each of the different solutions faces problems within a particular subset of the data. 
The ’3+1’ mass spectra are in disagreement with the short-baseline experiments and the ’2+2’ 
mass spectra conflict with either the atmospheric or the solar neutrino data. Therefore the 
four neutrino models, though not completely ruled out at present, seem highly disfavored. 

In the absence of a sterile neutrino, Yanagida and Murayama have recently suggested a 
another possibility to solve all of the known neutrino anomalies . The Yanagida-Murayama 
scheme preserves Lorentz invariance^ but involve CPT violation by invoking independent 
masses for neutrinos and antineutrinos. Hence there is a total of four independent Am^’s. 
Schematically we can represent the masses for the neutrino and antineutrino as in Fig.^ The 
solar neutrino problem only concerns the disappearance of z/g and the LSND experiment sees 
Vf. appearance in a beam. These experiments can be separately explained by Arug and 
AmLSND Yanagida-Murayama scheme. The atmospheric neutrinos data involve both 

and and allowing for CPT violation the two Am^’s are no longer constrained to be identical 
as also considered in Refs. IH] , ^ . It is therefore clear that all the data can be explained within 
this model. 

To invoke CPT violation is indeed a very drastic solution. Therefore it is important to 
discuss physical models for generating CPT violation in the Yanagida-Murayama scheme. In 
Ref. P] a dehnite model of CPT violation is introduced that can account for all available 
neutrino data. It was argued that CPT violation in the neutrino sector can be motivated from 
string theory via the extra dimensions. The right-handed neutrinos, like the graviton, are free 
to propagate in the bulk, whereas the Standard Model fields are constrained within a four 
dimensional brane. This gives rise to non-locality for the neutrinos and thereby generating 
CPT violation. This CPT violating scheme is also able to account for baryogenesis in a natural 
way 10. Furthermore non-commutative geometry can generate CPT violation [^, 0. 

In this paper we will study the atmospheric neutrino anomaly within a two family neu¬ 
trino scheme with CPT violation. The electron neutrinos are assumed not to oscillate on the 
atmospheric scale. We include the data from the K2K long baseline experiment, that further 
constrain the neutrino parameters. The atmospheric neutrino problem is by now well estab¬ 
lished and can be accounted for primarily by a two neutrino z/^ —> Ut oscillation 0. However, 


has been argued that the scheme also violate Lorentz invariance 
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sub-dominant oscillations are still possible and maybe even welcome [0, |TI|. Having different 
mass matrices for neutrinos and antineutrinos naturally gives different mixing matrices; IJi, 
for the neutrino sector and Ujf for the antineutrino sector. We will investigate whether the 
mixing parameters can be constrained by the atmospheric and K2K data. A most relevant 
parameter is the difference in mass squared difference for neutrinos and antineutrinos. Let us 
define the parameter e to describe the amount of CPT violation 




- I 

z>',atm I 


( 2 ) 


Using the latest data, we will show that e is only weakly constrained. 

Let us hnally mention that the LSND result, which has not yet been conhrmed, will be 
scrutinized by the Mini-BooNE experiment m- However, as has been noted before, unless 


this experiment is done also with antineutrinos the Yanagida-Murayama scheme cannot be 
ruled out. 


2 Analysis of the atmospheric data 


A number of experiments have measured the atmospheric neutrino fluxes. Here we will only 
consider the contained events of the Super-Kamiokande (SK) experiment |^. The justification 
for leaving out other data sets is the superior statistics of the SK data. Furthermore, the high 
energy upward through-going muon events [Wl are less affected by antineutrinos. For the 


average energy of 100 GeV of these events the flux ratio is about 1.5, thus decreasing 

the influence of the antineutrinos. Also the statistics is lower and we do not expect large 
effects from the inclusion of this sample. 

We use the following simple two-family survival probability relations for neutrinos and 
antineutrinos 


P„ 


Pv 


= 1 - sin^(26'j,) sin^ 


/ LAml' 


= 1 — sin^(26*i; 


. 2 f LAim^\ 

V ) 


( 3 ) 


( 4 ) 


We assume that the oscillation is into r-neutrinos, whereby the electron survival probability is 
taken to be one for both neutrinos and antineutrinos. As we only consider to Vr oscillation 
there are no matter effects. The pathlength of the neutrino, L, is calculated using an average 
production point in the atmosphere of 15km. E is the neutrino energy. 

The data are divided into sub-GeV and multi-GeV energy ranges and can be represented 
as the ratio, between the measured fluxes and the theoretical Monte Garlo prediction in 
the case of no oscillation. We define as 


10 

= E E E 

M,S a=e,ii i=l 


{Raj - Rj,r 


a" 


+ x} 


( 5 ) 


where aa,i are the statistical errors and M, S stand for the multi-GeV and sub-GeV data 
respectively and i denotes the zenith angle bin. For the details of the dehnition we re¬ 
fer to Ref.[|^, except that we here use a smearing of the sub-GeV events with an angle 


50°/ Eu/GeY [^|. The overall normalization of the neutrino fluxes is allowed to vary freely. 
Hence we minimize with respect to a, where the neutrino flux is given by $ = (1 -f The 

theoretically predicted neutrino flux is taken from [^, The y| term takes into account 
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Figure 2: The 68.3% (90%) CL regions for parameters Ami Am|- for the contained SK 
events and the total number of events in K2K (5 d.o.f). 


the error in the flux ratio. The SK Collaboration estimate the error to be 8 % in the 

sub-GeV range and 12% in the multi-GeV range. We renormalize the neutrino fluxes as 


|,s,Ar 




^S,M 

e 


(l+/3s,M/2)4f'" 


( 6 ) 


where the symbols S', M stands for sub-GeV and multi-GeV respectively and minimize the 
total function with respect to j3s and /3 m- The function is given by 




0.08 


+ 


0.12 


(7) 


For the best fit point obtained by the SK Gollaboration the values used are Ps = 6 % and 
/3 m = 12 %, implying that while scaling the electron ratios down one simultaneously scales the 
muon ratios up. 


We also include the recent data from the K2K long baseline experiment 


The beam 


is almost pure and we will neglect the small contamination of and z/g. We use the same 


method as in Ref, and only fit to total number of observed events. In total we have five 
parameters {Ami, 9,^, 6^, a) and 41 data points. 

The minimum is xlain = 33 at a = 1%, Ps = 8%, /3m = 10% and 


Ami = 2-5 X 10 “^ eV^ , Am| = 2.0 x 10 “^ eV^ , sm^{2e^) = 1.0 , sm'^{2eTj) = 1.0 ( 8 ) 


with 36 degrees of freedom. This is very close to the CPT conserving case. In Fig.^ we show 
the 68.27% and 90% conhdence levels as obtained by Ay^ < 5.9, 9.2, respectively, for hve 
degrees of freedom. At 90% G.L. the mass squared difference for neutrinos is constrained within 
4.5 X 10“^ eV^ — 5 X 10“^ eV^, while the antineutrino mass squared difference is only bounded 
from below (> 2 x 10“"^ eV^). At 90% G.L. the mixing angles are bounded, sin^(26'i,) > 0.8 
and sin^( 26 * 17 ) > 0.5, but maximal mixing is preferred for both angles. 
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Figure 3: Predicted ratios as a function of the zenith angle for mixing parameters; Ami ~ 
2 X 10“^ eV^, Am|. = 0.1 eV^ and maximal mixing, and the data with statistical errors. The 
upper curves are the ratios and the lower curves are the ratios. The triangles are the z/g 
experimental ratios and the circles are the experimental ratios. Note that the experimental 
ratios are plotted using <l> and not <F°. 


Multi-GeV 




The correlation between the lepton and the neutrino angle in the sub-GeV range is very 
weak and the data is smeared compared to the multi-GeV sample . The exact calculation 
method can therefore change the results slightly as the effect of a large Am|. is similar to a 
smearing. For smaller effective smearing the best £t point will move toward larger values of 
Am|.. The same considerations apply to the mixing angles which also effectively flatten the 
zenith angle curve. 

One should note that a point with values of Am|. ~ 0.1 eV^ and Ami — 2 x 10“^ is not 
signihcantly disfavored. The SK contained data become independent of Am| in this region 
as the oscillation probabilities are averaged to 1/2 for all pathlengths. K2K is obviously 
independent of Am|. as they measure neutrinos. The predicted ratios of this mass spectra 
are shown in Fig.^ along with the data points. The predicted ratio of around 0.85 for the 
downward going multi-GeV muon neutrinos can be easily understood. In this energy range 
the flux of neutrinos is roughly the same as the flux of antineutrinos. But the antineutrino 
cross section is less than half that of neutrinos. The sin^(LAm^/4F^) is averaged to one half 
for antineutrinos and to one for neutrinos and the ratio is estimated to 


i? I ~ 


$ P 


,cr„ 


<^_ P— — (J— 


i'/J, 


~ 0.85 , 


(9) 


where a is the cross section. The reason that this mass spectra is not strongly disfavored is 
because it agrees very well with the well known double ratio. The measured value is 


R = 


(f/g 


DATA 


(h/e) 


MC 


= 0.675iS;S| ± 0.080 


( 10 ) 


The prediction double ratio for Ami — 2 x 10 ^ eV^ and 
Am^ ~ 0.1 eV^ are R ~ 0.68, whereas for both Am^’s around 3 x 10“^ eV^ we get R ~ 0.75. 


in the multi-GeV range [IT8 
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In both cases we have assumed maximal mixing. For sub-GeV events the double ratio is 
0.638 ± 0.017 ± 0.050. In other word the advantage of the CPT violating mass spectra is that 
it for a flux normalization that diminishes the excess of in particular sub-GeV Ue events also 
agrees well with the muon ratios. In fact in the case that the overall normalization is varied 
freely, but the /r/e ratio is kept hxed, the best £t point is for Am| Ami- Although such 
a mass spectra does not fit the up-down asymmetries that well as seen from Fig.^. It must 
be remembered that the measured values quoted in Eq.(p!0[) is dependent on the theoretical 
predictions for the z/^/z/g flux ratio. The double ratio basically describes the average gap 
between the and the z/^ ratios and this gap is somewhat too large to be fitted very well by a 
CPT conserving two family z/^ —oscillation. However if the theoretically predicted v^. 
flux ratio is decreased by 6-12% this provides a very good fit to the data. Furthermore the low 
value of the double ratio could be due to an excess of Vf. events which is not accounted for by 
the two family scheme. The value for the point in Fig.|^ is 39 and therefore barely outside 
the la region. When using 2 d.o.f. as in Ref. p6| this would be a 2a exclusion. We remark 


that there is a local maximum of the function for values of Am|- between 10“^ — 10“^ eV^. 
The main differences between our results and those in Ref. [Q are due to the fact that we use 
the theoretical predicted fluxes, whereas the normalization in each type of events are varied 
freely in Ref. ||26|| . 

For almost all mixing parameters a rise of the theoretical fluxes are needed. New calcula¬ 
tions suggest instead a lower flux, by about 10% [^, mainly due to the primary flux being 
lower than obtained in earlier measurements. The new normalization results in a large excess 
of z/g events which seems to be difficult to obtain theoretically. In the conventional CPT con¬ 
serving case the SK Gollaboration obtains a = 20% for the best fit point using the new ffux 
predictions |^|, which roughly amount to putting the normalization back to the old value. 


Therefore we find the error of the overall normalization is large and a rise can not be excluded. 

In the CPT violating scheme the electron ratios can be considerably away from one. A 
relatively large LSND angle [0 as well as the solar mass squared difference [^] can inffuence 
these ratios. Though the LSND angle is constrained to be small by the BUGEY results. Also 
the two angles 6*g^ and 9^^, describing the oscillation of z/g driven by the atmospheric mass 
squared difference, can have effects. In particular 9’^^ 
and Palo Verde 


is not constrained by the GHOOZ 
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results as these experiments are measuring anti neutrinos. Moreover 
can be large if AmZ is below the GHOOZ sensitivity of 10“^ eV^, which is not ruled out by 
the present data (see Fig.^). The influence of some of these extra mixing parameters has been 
studied in Ref. |P2||, where however the systematic errors in the SK data have been ignored. 

As we have shown the limits on CPT violation in the neutrino sector are rather weak 
at present. There are nevertheless good prospects for a much stronger bound in the near 
future. The results from the KamLAND experiment [^] could likely disprove the Yanagida- 
Murayama scheme. The experiment will test the currently favored large mixing angle MSW 
(LMA) solution to the solar neutrino problem, by detecting anti electron neutrinos from nearby 
nuclear reactors. For the most favored region of the LMA solution (Am^ < 2 x 10“"^ eV^) the 
detected energy spectrum at KamLAND will quite precisely determine the value of the mass 
squared difference and this signal would rule out the Yanagida-Murayama scheme. However 
for Am^ > 2 X eV^ the oscillations are averaged out and KamLAND can only put a 
lower bound on the anti neutrino mass squared difference In this case the situation 

becomes more problematic as there are two different possible explanations. A large value 
of AtUq is not ruled out, though disfavored by the present data. The signal could also be 
explained within the Yanagida-Murayama scheme by having a small along with a 

large 0g^. Borexino detecting solar neutrinos will not be able to pin down the true solution. 
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Hence, if an averaged oscillation with a large angle is observed, one would most likely have to 
wait for the results of the MiniBooNE experiment. In the case that KamLAND does conhrm 
the LMA solution to the solar neutrino problem, much better limits on CPT violation could 
be obtained as discussed in Ref.|^. If Kamland does not observe a suppression of the anti 


neutrino flux there are different possibilities to test CPT violation as discussed in Refs. [26, 32 


Indirect limits can also be obtained by studying how radiative corrections communicate the 
large amount of CPT violation in the neutrino sector to the charged lepton sector 


36 


In conclusion we have analyzed the Super-Kamiokande contained events and the K2K 
data in a CPT violating two neutrino z/^, framework. The best fit area is close to the CPT 
conserving case. However at present the CPT violation parameter e, defined in Eq.(|^), cannot 
be usefully constrained. 
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